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SUMMARY 


A group of problems involving time -dependent discharge of 
inviscid liquid from a reservoir is examined. Two geometries are 
treated in detail: the slender conical reservoir, and the cylin- 
drical reservoir. 


The differential equation governing free-surface motion 
during discharge is non-linear and cannot be integrated exactly 
except for the special geometry of a cylindrical reservoir. How- 
ever it is shown that it is possible, in the two specific cases 
treated, to neglect the free-surface acceleration. This is called 
the approximate unsteady solution. This is compared to a numeri- 
cal solution of the exact equation, the quasi-steady solution, and 
experiment. It is highly probable that free-surface acceleration 
plays no role in the flow regardless of the reservoir geometry. 
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For the cylindrical reservoir it turns out that only with an 
exit nozzle is it possible to make meaningful comparison with ex- 
periment, and even then special precautions must be taken to 
avoid free-surface collapse. The correction to the discharge time 
even for short nozzle lengths may be a significant fraction of the 
total. Experiments indicate that when the length/diameter ratio of 
the nozzle is large, friction in the nozzle is important. When the 
ratio is small the inlet shape must be carefully designed in order 
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to avoid flow separation of the jet from the nozzle walls. 

Other theoretical work is discussed briefly. In particular 
the work of Stary is considered who employs momentum and kinetic 
energy correction factors to account for the lack of uniformity 
of the flow across any cross-section normal to the axis. It is 
shown for a particular geometry that the introduction of these 
correction factors is superfluous and it is concluded that their 
significance even for steady flows should be reexamined „ 
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F, G 
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cross-sectional area 

inverse of velocity coefficient 

diameter 

defined in eq, (19) 
gravitational constant 
Stary's correction factors 
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indices or exponents 
pressure 
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vertical velocity component, 
dummy variable 

free-surface coordinate, dimensionless 
vertical coordinate, dimensional 
A2 / A^ 

defined by eq. (15) 

z 2^ z l f° r c yl indrical reservoir 

free-surface coordinate, dimensional 

z 2^ z l ^ or COT1 i ca - reservoir 

fluid density 

dimensionless time 
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INTRODUCTION 


The subject of unsteady liquid flow through ducts has net 
been greatly developed. For example, in the recent ’’Handbook of 
Fluid Dynamics", reference 1, no mention is made of the subject. 

It is possible to deal with flows with negligible viscous effects 
by standard methods of potential theory involving superposition 
of singularities. Roudebush and Pinkel, reference 2, treat the 
unsteady flow out of a cylindrical tank by this method but it is 
too complicated to have much promise except in special cases. 

One reads the standard works in vain for an adequate treat- 
ment of the subject. In their famous work Prandtl and Tietjens, 
reference 3, approach the subject by means of the unsteady 
Bernoulli equation. They point out (correctly) for one situation 
that the governing differential equation is of second order but 
no details are given. 

Kozeny, reference 4, in a work available only in German, 
uevelops the theme outlined by Prandtl. In a short note Stary, 
reference 5, applies the development of Kozeny to the problem of a 
cylindrical tank with a horizontal bottom which is provided with a 
discharge hole. Neglecting viscosity he obtains a solution for 
the discharge time in terms of the area ratio, the original free- 
surface height, and certain momentum and kinetic energy correction 
factors which have to be determined empirically. Bird, Stewart and 
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Lightfoot, reference 6, page 239, pose as an exercise, the same 
problem considered by Stary; their expression (taking into account 
what appears to be a misprint) for the discharge time agrees with 
that of Stary when his correction factors are put equal to unity. 
We shall also give the solution to this problem by an alternative 
approach. 

Bird et al, also consider, page 226, the unsteady discharge 
from a conical tank, which is a special case of a class of flows 
considered in the present work. By neglecting the kinetic energy 
they obtain a result which is equivalent to a quasi-steady flow 
whose exit velocity is given by the Torricelli value based on the 
instantaneous free-surface height. 

Figure 1 illustrates the conical reservoir. Assuming a 
fluid of negligible viscosity under action of gravity. we ask: 

If at time t = 0 the valve at the exit is opened, what is the 
discharge time t<j required to empty the reserx^oir? 

The quasi-steady solution. According to Torricelli's law 
the exit velocity at station (T) is 

= - '/tyX-T) ~ - fsfi, ( 1 ) 

if Zq ^ 3* • From the continuity equation for an incompressible 







k - i 
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fluid the velocity of the free-surface £ =j£(t) is w(j£) = d3*/dt 


2 , Jl 


w l z l » an< ^ thus the discharge time is 


£ “ ^ 






W-?J 


1 /* ) 


( 2 ) 


It is convenient to work with dimensionless quantities. Putting 
S (g/ z 2 )^t , )\ = z 2 /z 1} then 


£j ~~ A ^/< 


( 3 ) 


It is noted that according to quasi-steady theory the discharge 
. speed depends only on the free-surface height. Thus the discharge 
velocity theoretically becomes zero as the free surface reaches the 
exit. On the other hand in a non-steady incompressible flow the 
entire body of fluid must be accelerated from an initial velocity 
of zero. In a non-dissipative flow the free-surface velocity ought 
to increase in time and be a maximum at the exit. Therefore it is 
not obvious in advance whether the quasi-steady or ncn-steady 
theories ought to predict the lesser discharge times. 


Comparison of the solution of the unsteady flow equations with 
the quasi-steady solution will provide a convenient measure of the 
importance of non-steady phenomena. It will be shown that in some 
~ases u£e of the quasi-steady result can lead to significant error. 

THE UNSTEADY EQUATIONS 

-We assume that the flow can be treated as one-dimensional. 


There is no particular difficulty, initially, in generalizing tc 
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| include flow geometries involving ducts with area distributions 

following the law A = kz 11 . For axially -symmetric ducts n = 2 

t 

•- designates a conical boundary, whereas r. ^ 2 correspond to convex 

* and concave interiors respectively (see Figure 2). Furthermore, 

it is simple to include in the initial analysis the effect of a 
pressurized reservoir such that the free-surface pressure remains 
constant 5 

For unsteady incompressible flow the continuity equation is 

j urCzJJ/Jte) - 1 )/?&,) - (4) 

Thus, specification of the velocity at any station automatically 
\ specifies the entire flow field, 

I 

The dynamical equation is 

I - -J. - s. (5) 

\- /O 22 * ' 

: Specifying a constant discharge pressure p^, and constant 

c 

t 

free-surface pressure p 2 , the boundary conditions are 

l _ ' 

• 7° (?,,*) = /? j ' /°[W ] = . (6) 

£ 

s v 

! Initial conditions are 

l ' . 

i ~ 

f(o) -? 2 , «rfr s> 0) = *f[tfo)l= 0 . ( 7 ) 
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Equation (5) can be rearranged and integrated from the dis- 
charge station to the free surface at an arbitrary time 

J ^ * 

f 2 £ -f- f (■ -2- — + ex') 'Zyz- — n 

2>2 T / 2 ** ( ^ r 


2=2, 


t J 3k (p+ f +/*) >* = 0 


2 = 2, 


r 

+ **(*) i-?v-z>)=°- w 

* 


From (4) 


4) - 


4” - xr7*”. 


where 'f =. d//dt* Therefore 


•* *) - y? ' / . / 7? 

Jp = (S3 + r?3 3 )/2 2 


(10) 


3 

s' h? 

J -* 

*/ 


•• * *r- M-0 - M 1 / 

(33 + ??3 3 )[* - f J /'*")_ (11) 


for 'nj^l. 


Combining (8) and (11), we obtain the following differential 
equation for the free-surfade motion: ' - 


-/ .. Q-f -j r . £ fZ-f, f 

(7,-,) (*3+»s )[(■£) J/ + £_' +f(s-z,)=0' i 


( 12 ) 
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Now we define 


X 2 


s zjz, , jf * , T 3 (j/t t ) f y 


(13) 


and equation (12) becomes 


S M V ’‘.fe/ otf”-'] r^r-'p%+ p = < 14 > 


with 


/? = (fk-ffr/fjU 


(15) 


For the. conical duct of Figure 1, n = 2, and equation (14) reduces 


3 <? 


+ i [ 3 - W - ty)]} +i'+rfQt-i)-o. (i6) 

For n = 1 the integration of (10) introduces a logarithmic 
term. The final differential equation for this case is 

f tyyjf + i[i - +/ 3 _+ ty-'jfi = o. a?) 

Note that (16) and (17) are independent of the constant k; i.e. 
in the case n = 2, the flow does not depend on the actual cone 
angle, and similarly for other shapes. 


The governing initial conditions in all cases are 


= / , j(o)=o. 


(18) 


The solution terminates at y , when the free surface reaches 



the exit station. 


Equations (14) through (17) a^e surprisingly formidable 
for what is apparently a simple flow situation, i.e., a flow 
which is one -dimensional , incompressible and inviscid. A general 
solution to (14) has not been obtained. Therefore we focus 
attention on the conical duct of Figure 1, for which n = 2, which 
is governed by equation (16). 

INTEGRATION OF (16) 

A variety of attempts to integrate (16) were not entirely 
successful. It is useful to mention some of these briefly be- 
cause they point up the mathematical difficulties. Furthermore, 
it was in the attempt of trying to evaluate the utility of one 
of these approaches that the clue- to a successful approximate 
solution was obtained. 

Quadrature . Equation (16) can be integrated by quadrature. 
Putting 

ry) = y '' Ls - tyf - tyf - vJ ] 

► 

SCf) = 

0 0 n 

and noting th't y a d(%y z )/dy, then (16) becomes 


( 20 ) 
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Multiplying by the integrating factor exp j* F(y)dy, then 


first integral gives 


JLu* = 
2 


<f>[- f <■ 


and finally 


z, = 


/ X 


Evaluation of (22) depends on the integration indicated 
in (21) which has not been achieved. A numerical integration 
looks possible using a digital computer although great care 
must be employed because, in practical cases where values of 
2 C 100 might be encountered, extreme variations of the 
integrands can be expected. For this reason it was decided to 
investigate other possibilities first. 


Series solutions . Frobenius 1 method in which y(X) is repre- 
sented as an infinite series in % fails because the number of 
terms required for adequate accuracy (values of as large as 
2000 are encountered) would be very large. And then there is no 
guarantee that the series would be convergent. 


Asymptotic expansion in powers of a small parameter . An 
attempt was made, for ^9 = 0, n * 2,_ to put 
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substitute in (16), and solve for the resulting sequence of 
differential equations in y m . Unfortunately these differen- 
tial were incompatible with each other and the method failed. 
It seems likely that some variation of this technique ought to 
work but it remains to be discovered. 


THE METHOD OF SUCCESSIVE APPROXIMATIONS. 

Equation (20) is particularly suitable in form for employ- 
ing a well known approximate method for solving non-linear 
differential equations. One of the lower derivatives is re- 
placed initially by an expression which allows the resulting 
equation to be integrated for the term which was originally 
replaced. This expression is then put into the original differ- 
ential equation and, hopefully, the succeeding differential 
equation can be integrated. The process, if convergent, can be 
repeated until the actual solution is approached to any accuracy 
desired. 


Thus equation (20) can be rewritten and approximated as 


follows: 


tt 


/ 

where the subscript m denotes the order of the approximation. 
For n a 2 and ^ s 0 this becomes 


( 23 ) 








- 10 - 
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-\M - f £ C , f(s) cij . 

(y - • ton- ! 

/ 


(24) 


In principle we are indifferent to the choice of y Q as 
long as the process converges. Expecting that the exit velocity 
rapidly builds up it seems reasonable to try the Torricelli 
value, based on the initial free-surface height, as the zeroth 
approximation,, However, substitution of this value, which is 
y Q - y /2 A » in (24) leads to imaginary values for . This 
Indicates that the procedure is very sensitive to the choice of 
the zeroth approximation. 


On the other hand, choosing y 0 = 0 leads to 



125) 


Substitution of (25) to obtain the next approximation yields 



~ + ' v 7 + <2«) 

Unfortunately equation (26) also predicts imaginary values 

* 

for j 2 ancl the method again fails. Furthermore whe complexity 
of (26) indicates, even if a satisfactory zeroth approximation 
could be found, that if convergence were slow, requiring several 
iterations, the method might still be unwieldy. Furthermore, if 



11 


an expression for the velocity could be obtained with adequate 
accuracy, evaluation of the discharge time would still require 
an integration in which the last approximation appears as a 
radical in the denominator of the integrand. Integration, ex- 
cept by a numerical procedure, of such a result appears unlikely. 

This infelicitous behavior is not without some redeeming 
consequences, however. It is seen from (16), for - 0, that 
initially, since y = 0, the second order term is dominant. How- 
ever, as the free-surface speed builds up, the magnitude of y 
must decrease rapidly from its initial value of A because the 
coefficient of y^ is initially of order Since y cannot 

change sign this suggests that the second order term may in fact 
be negligible throughout most of the discharge. Before investi- 
gating this point a criterion for the time to "start' 1 an unsteady 
flow is propounded- 

THE STARTING TIME 

There is no unambiguous way to define a starting time for 
an unsteady flow. At the exit the velocity is initially zero. 
Intuitively we expect that the velocity rapidly builds up until 
it reaches or surpasses the Torricelli value. When this condition 
is reached the free-surface velocity is y = -Jl A'^for n = 2 
and A - 0, We define the starting time ^ as the time to reach 


this value. 
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An approximate lower bound on can be obtained as follows. 
Expecting that the free-surface build-up is attained rapidly we 
put y«l, y ' 2 Z 0 in (16), and obtain y~ which is integrated 

to give yc* ' This is equated to the Torricelli value for 

y = 1 and solved for the starting time which is 


5 


Si\ 


-/ 


(2 7) 


Values of 7T S , according to (27), turn out to be of order 


-3 


^ times the discharge time for the same configuration, which 
is an insignificant fraction of the total. Similar expressions 
can be derived for arbitrary n and ^3 , 

A N APPROXIMATE RESULT FOR THE DISCHARGE TIME 

For n s 2 andy3 * 0 the free-^urface acceleration is 
initially, y = As y^ increases, y decreases rapidly in 

magnitude and remains small for most of the discharge run. As 
the free surface approaches the exit, however, the volume :>f 

liquid decreases at a rapidly increasing rate until for y * , 

» • 

y = -1, which corresponds to the acceleration of gravity., 

This suggests that in (16) it is possible to ignore the 
acceleration of the free surface, reducing the problem to the 
integration of the following first-order equation: 


4 


[3- ty) 3 - ty) - \#]i 


-f~ \ — o , 


(28) 
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It was first pointed out by Prandtl, whose discussion is 

reproduced by Schlichting, ref. 9, page 63 5 that it is not 

generally possible to neglect the highest order derivative in a 

differential equation because this means that one of the initial 

conditions cannot be satisfied. In the present case, however, 

the approximation leads to reasonably accurate discharge times. 

-2 

It is conjectured that from (28) although y(0) » , which 

is the Torricelli value, it is nevertheless sufficiently close 
to the exact requirement y(0) = 0 that the resulting error for 
the discharge time is negligible. Equivalently, use of (28) is 
tantamount to neglecting the time to start the flow as given by 
(27). 

On the other hand (28) breaks down when the free surface 
reaches the exit, yielding an infinite velocity at that point. 
Since the time for the free surface to cover the last small 
fraction of its total displacement is small in any case, the 
singular behavior at the exit has no great effect on the dis- 
charge time. 

Solving (28) for the negative root of y and integrating 

gives 

Ti = /vi / [ty) + W + V ~ 3 1 ‘‘ti < 29 > 
d y 


or, equivalently. 
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Since the second term in brackets is small compared to the first, 
except near the exit, we ignore it to obtain the explicit result 


% -- \ z ji/s- -t - (31) 

\ 

keeping only terms in / to the first power or greater. 

The first term of (31) is the sane as the quasi-steady ex- 
pression of equation (3). Consequently the solution for the full 
unsteady flow problem is a second-order correction to the quasi- 
steady solution. This is borne out by Figure 3 in which are 
plotted equation (3), equation* (31), a numerical integration of 
(29), and experimental values -- all for n = 2 -- and for n = 1, 

a numerical integration of equation (17) for ^3 - 0, neglecting 
the term in y. 


During the preparation of this paper equation (16) was pro- 
grammed for the Philco "Transac n digital computer and several 
runs were made. For the case \ = 40, /S = 0 the result 
differed only by .one part in 900 from that of equation (31). 
For ^ = 10, 10 the digital computation yielded ^ s 7.24, 

which is about 1% greater than that of a numerical integration 
of (16) with y s 0, It is amusing to note that in the digital 
computation the function y oscillated irregularly in sign all 
the; way to the exit although its magnitude averaged only about 
10“ . This behavior caused only a few wiggles in the velocity 
term y, and can be taken as another indication that y can be 
safely neglected. 
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There is no significant difference between the values given 
by (29) and (30) except for X ^ 5. For ^ - 10 the value given by 
(31) is about 8% higher than the quasi-steady value, for A = 40 
this drops to 2% and to 0,8% at An 100. 

Experimental data were ''btained on equipment, described in 
the Appendix, which permitted use of any one of three different 
orifices. By varying also the liquid height, the range of values 
6 £ A— 50 was obtained. Experimental values for were slightly 
higher than predicted by (31) varying from 3% in most cases to 10% 
foi: a few. ” It is reasonable to expect that a part of the differ- 
ence between theory and experiment is due to the neglected action 
of viscosity. 


THE CYLINDRICAL RESERVOIR 


Theory , The flow from a cylindrical reservoir with a small 
nozzle at the exit has been posed as an exercise by Bird, Stewart 
and Lightfoot to be treated by one-dimensional theory. Their 
result for the discharge time is given in the form of an unevaluated 
integral. Stary in reference 4 tackles the same problem using the 
theory developed by Kozeny, in which the unsteady Bernoulli equation 
for a streamline in a general potential flow is integrated first 
transversely and then streajnwise, yielding a differential equation 
for the free-surface motion which is integrated to give an expression 
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for the discharge time. His expression involves certain empiri- 
cal factors intended to correct for the non-uniform distribution 
of kinetic energy and momentum flux when integrating over a cross- 
section. Surprisingly, no experimental results appear to be extant 
for this flow, so that it is not possible to evaluate either treat- 
ment . 

It is desirable to compare the one-dimensional treatment with 
that of Stary to evaluate the importance of his correction factors. 
A schematic is shown in Figure 4a We commence with a concave in- 
terior connecting the fixed cross-section area A 2 , located at 
z = Z 2 , with an outlet area A^ at z = 0. For the initial area 
distribution we choose 

/9 = /?/// (°'-/)(s/? s )’] } (32) 


where = A 2 /A 1 . Omitting consideration of a pressurized reservoir, 
we repeat essentially the same procedure by which (14) was derived 
and are led to the following equation for arbitrary n: 
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Taking the limit as n-*>0 produces the desired cylindrical 
contour, for which the governing differential equation is 


y-f + i (/-<=<■) if- + p - o, 

with the same initial conditions given by (18). Equation (34) 


( 34 ) 


can be integrated as it stands in terms of gamma functions. We 
define s = at -1 and then, employing the transformation y = u , 
we obtain a form whose integral is given by Dwight, reference 9, 
formula No. 857.1: 


*? - 

O u J 


& r[sjh>] 

= [ — 

for s > 1; for s s 1 


(35) 


Kl = 


(36) 


An expression similar to (35) is obtainable for s < 1, but this 
case is of little practical interest except to note that in a free 
fall which is the limit as <>£-*1, -rjl* We have verified that 
the expression of Bird, Stewart and Lightfoot, taking into account 
a misprint, produces these same results assuming a negligible nozzle 
length. *. 




»*<.> „ -- . TfW :/ >**««• < -■ • r>r n r - . i_V. 





The approximate unsteady solution, neglecting the second- 
order derivative in (34), is 

'Q = < 37 > 

which for 1 reduces to the approximation 

- /F oC. (38) 

Equation (38) is also the result obtained from quasi-steady theory. 
Obviously, the approximate theory of (37) breaks down for area 
ratios near unity since it predicts £^-*0 when o^-*l. For area 
ratios c/ ^ 10, equations (3 r >), (37) and (38) are indistinguishable 
for practical purposes. For: example when oC - 10 exact theory 
gives = 14.10, while the approximate and quasi-steady values 

are 14.07 and 14.14 respectively. We conclude that for all area 
ratios of practical interest, unsteady effects are negligible. 


The quasi-steady theory starts to diverge somewhat from the 
exact theory below oC » 5 but curiously, for free fall, when ^ - 1, 
they both yield = J~2 , This is all the more remarkable when 

we consider that the free-surface velocity according to quasi-steady 


* After completion of this analysis it was brought to our attention 
that Kaufman, in a book reference 8, recently translated from the 
German, has also treated the problem of reservoir discharge. For a 
general area distribution he obtains a quadrature similar to equation 
(22), and then for a cylindrical container he puts y « 0, and obtains 
eq. (38) modified by a velocity-coefficient correction factor.. 



theory is wrong at every instant except one, and that it varies in 
the wrong sense with respect to time. The reason it yields the 
same discharge time is that, if in the free-surface velocity func- 
tion, according to quasi-steady theory, we substitute ^ - 2 T for 
'fc we obtain the correct history in a free fall. Integration of 
either function of course must yield the same discharge time.. 

Discussion of Stary 's solution . By neglecting the free- 
surface acceleration Stary obtains [ eq . (4) of his note rewritten 
in the present notationj ^ = £2(kio4“* - k 2 )j^si k^ \TloL where 
k^, k 2 > k .3 are correction factors, and where 1.04 £'^3 1.08. He 

does not explain how to choose the appropriate value for k 3 , Stary 

„ _ 7 

terms this the quasi-steady solution. If we put k 3 = 1 then his 
result reduces to (38) , which we have called the approximate un- 
steady solution. Of course for the cylindrical reservoir there is 
little distinction between the two. However for the conical reser- 
voir the quasi-steady solution underestimates the actual discharge 
time which, on the other hand, is given quite accurately by the 
approximate theory. In other words neglecting the free-surface 
acceleration is not equivalent to the quasi-steady approximation. 

Stary also obtains a solution of the complete equation which 
gives the. ratio of the exact discharge time to the approximate 
discharge time. The ratio is less than 1.02 for 0 ^ 2 10 and approaches 
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unity as the area ratio becomes large. If the correction factors 
are set equal to unity then his solution reduces to equation (35). 
Thus except for the correction factors the unsteady Bernoulli 
equation and the one -dimensional theory are equivalent. 

The sharp-edged orifice . So far we have said nothing regard- 
ing flow conditions at the exit. For the conical reservoir, as 
long as the cone is slender, the contraction effect of the jet is 
small as our experiments show. For the cylindrical reservoir how- 
ever, with a sharp orifice at the exit -- the problem considered by 
Sfary -- the jet leaves with a significant radial component of velo- 
city inward which causes a radical contraction of the jet area 
downstream of the exit. The effect on the discharge time is enor- 
mous as shown by Figure 6 in which Stary's approximate solution, 
with kg = 1.08, is compared with experiment for an initial free- 
surface height of Zg = 11.25 in. "'"he discrepancies are too great 
to be attributable to lack of precision in the correction factor 
k^. Instead, following the usual procedure, it might be better 
to put ^ = C y~2ot where C is the inverse of a velocity coeffici-' 
ent, and is evaluated from the experimental data. Whether or not 
it is valid to employ a non-time-dependent velocity coefficient 
in an unsteady flow is not known. 

Whether or not to treat an unsteady flow situation by quasi- 
steady theory depends on factors which a.re not easy to pinpoint. 



Sabersky and Acosta, in a recently published text book reference 
10, treat the flow out of a cylindrical reservoir by quasi-steady 
theory and conclude rightly that unsteady effects are negligible 
for area ratios 2 10. On the other hand their analysis is 
equally applicable without alteration to a reservoir of arbitrary 
geometry, a conical reservoir for example. As we have shown, even 
for area ratios up to AO, quasi-steady theory for the conical 
reservoir predicts discharge times from two to eight percent less 
than the approximate theory which itself is less than experiment. 

Thus it is our opinion that their method for determining a crite- 
rion for neglecting unsteady effects is faulty and should be used 
only with caution. 

Theoretical correction for the nozzle length . Rather than 
introduce the additional complexity of a velocity coefficient we 
have' elected to deal with the configuration of Figure 4b where a 
short transition section and straight nozzle enable the jet to be 
emitted without any radial component of velocity at the exit. The 
velocity coefficient of a well-shaped, smooth nozzle can be ex- 
pected to be near unity except when the ler.^ h/diameter ratio be- 
comes large. In this event viscous effects may become non-negligible 
as can be visualized in the limit of decreasing diameter as the 
nozzle becomes a capillary tube. 



22 


For a short nozzle-plus -trans it ion-section length the asso- 
ciated volume of liquid is negligible with respect to the whole. 
Therefore for analytical purposes the flow terminates effectively 
when z = z^. By subtracting the dimensional discharge times for 
two flows of original heights Z£ and zj_, and then non-dimens ionaliz- 
ing using as the reference length, we obtain the discharge 
time, denoted ^ , for the nozzle configuration 

I- (/- £ ,& ), (39) 

where 6 = Equation (39) is based on quasi-steady theory 

but it is easy to show that it is also valid for exact theory to 
a close approximation, for values of 6 not too near unity. 

EXPERIMENTAL RESULTS FOR THE CYLINDRICAL RESERV OIR 

Although measurement of the discharge time for a cylindrical 
reservoir appears almost trivial to carry out there are several 
surprises for the unsuspecting. It is possible that failure to 
cope with these unforeseen effects accounts for the dearth of 
published data for this flow. 

Effect of nozzle shape . Early measurements with a straight 
exit preceded by a transition section which was simply a 90° 
circular arc section of 3/8 in. radius, resulted in discharge 
times substantially greater than predicted by theory, presumably 
due to flow separation from the nozzle wall. We then went to the 
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nozzle s hown in Figure 4b in which the inlet was shaped to contour 
A of Figure 5, This contour, which is not based on a theoretical 
design, was employed for all of the runs with the exception of a 
few at a = 3 in. We shall return to this point shortly. Details 
of the apparatus are given in the Appendix. 

Free-surface collapse . Another phenomenon occurs for the 
smaller area ratios (for d - 1 in. and above). As the free-surface 
height decreases, a point is reached where the center is visually 
lower than that of the liquid adjacent to the walls. This is a 
consequence of the pressure gradient normal to the streamlines, 
near the transition section, as the fluid makes a turn more or less 
parallel to the nozzle boundary. When conditions are favorable the 
free surface, which remains at atmospheric pressure, pops through 
the nozzle forming a long finite cavity inside the jet. This des- 
troys the almost one -dimens ional nature of the flow, decreases the 
effective jet cross-sectional area, and consequently increases the 
discharge time drastically. 

To eliminate this interesting but unwanted effect we employ a 
thin balsa wood disc, a little smaller in diameter than the reser- 
voir, and which floats on the free surface. Surface tension "seals" 
the free surface and eliminates the popping through entirely. 

Although this phenomenon does not occur for the smaller holes 
(d-^ - 0.25 and 0.5 in.), use of the disc is helpful there also. 
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As the free surface approaches the terminal level zj , surface ten- 
sion and friction become more important. The result is that, 
without the disc, instead of the flow terminating in a sharp cutoff 
it tails off gradually into a dripping process. With the disc 
however, there is a distinct transition from a continuous jet to 
the discrete droplet stage, which was taken to signal the end of 
the flow,. We have not heard of this device being exploited else- 
where but if there are industrial processes where cavity formation 
is a problem (such as in the draining of a large tank of molasses) 
it may have some potential usefulness. 

Turbulence in the je t. Even with the improved nozzle the flow 
in each case is initially turbulent on tne jet boundary. This is 
under ndable since the Reynolds' number, based on the initial 
quasi-steady velocity, and for water at 75°F, is 24,000 for the 
one-quarter inch hole, which value exceeds the minimum critical 
Reynolds number 2300 and is great enough almost to guarantee the 
appearance of turbulence. To the eye tne scale of turbulence in- 
creases as the exit diameter is increased. As the free-surface 
level falls the outflow Reynolds number decreases. Eventually a 
point is reached where the jet suddenly becomes laminar. The 
smaller the hole the sooner this occurs. 




Ref. 7, page 35 > 
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Experimental results. Experiment and theory are compared in 

Figure 7 , Runs were made for an overall nozzle-plus-transition 

length of z ^ = 1.25 in. and a height of Z 2 = 25.15 in., which corres- 

% 

pond to 6 = v.,0497 and 1-6 n 0.778. Each experimental point is 
the average of two runs. In no case does the deviation of any run 
exceed one percent of its corresponding average. Because a log-log 
plot tends to obscure the magnitude of error, a summary of our re- 
sults appears as Table I. 


TABLE I 

Experiment and Theory for the Cylindrical Reservoir 



z 2 

(in.) 

— 

6. 


*5 

Eq. (39) 

(sec.) 

(exp.) 

Dev. 

m 

Contour 

■ 3080 

1 

25 

.15 

0.04 

►97 

1 

0.25 

3384 

939.1 

3679 

8.7 



i 

| 770.1 




1 

0.50 

846.1 

227.8 

892.5 

5.5 


i 

i 

I 

! 192.5 

} 

t 





1.00 

211,5 

54.6 

213.9 

1.1 



| 85.56 





1.50 

94.00 

24.6 

96.38 

2.5 



i 

| 48.13 






52.87 

13.63 

53.40 

1.1 



1 

|. 30.80 






33.84 

9.25 

36.24 

7.1 

' 


21.39 






23.50 

7.90 

30.95 

31.7 

i 

\ 

21.39 


r 




23.50 

6,60 

25.86 

10.0 

B 

21,39 

25 

.15 

0.0^ 

►97 

3.00 

23.50 

6.60 

25.86 

10.0 

C 

3080 

13 

.25 

0.0943 

0.25 

3018 

611.3 

3309 


A 

192.5 

13 

.25 

0,0943 

1.00 

188.6 

35.1 

189.5 



0.6 

A 





Z 1 

= 1.25 

in. d 

2 * 13 • 88 in . 
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There is some disagreement between theory and experiment on 
both ends of the curve. For large values of the area ratio (small 
diameter holes) it was foreseen that viscous effects would be great- 
est. A crude idea of the magnitude of the viscous contribution can 
be obtained as follows: for an initial Reynolds number of 24,000, 
and a nozzle length of 0. 75 in. , theory of turbulent flow through a 
smooth pipe predicts a head loss of 1.9 in., this from an initial 
head of about 25 in. For the three-inch hole the corresponding 
Reynolds number is 290,000 and the head loss 0.09 in, A theory for 
variable Reynolds number is not developed, but Kaufman, reference 8 
page 123, gives an approximate scheme to hahdle it. 

As the initial Reynolds number becomes larger the deviation of 
experiment from theory becomes almost negligible and continues so 
over most of the central region of the curve. Thus we conclude that 
viscous effects are of importance only for the smaller nozzles , more 
specifically when the length/diameter ratio exceeds unity. We also 
conclude that there is no need to incorporate the kinetic energy and 
momentum correction factors of Kozeny and Stary in the analysis. 

Only by putting = 1 would Stary's solution agree closely with ex- 
periment for values of the area ratio near oL — 200. Furthermore, 
as we have indicated, for large values of cC the deviation is due 
to viscous effects and not the failure of the one -dimensional theory 
to account for variations of the kinetic energy and momentum flux in 
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the transverse direction. 

For = 3 in. , the largest diameter nozzle tested, uncorking 
the hole results in a rather startling gush of water. The theoreti- 
cal running time is only about six seconds. For contour A the 
actual discharge time was 30% above theory. For this configuration 
the jet was initially turbulent but turned laminar much sooner than 
for any other diameter. It is probable that the flow separated from 
the transition section wall although no attempt was made to verify 
this . 

Instead, another run was made using contour B, which has a 
slightly more, gradual transition terminating in a one-half inch long 
straight nozzle. This brought the experimental time to only 10% 
above theory. A third run, using contour C, resulted in no addition- 
al reduction in the discharge time. Since the theory must be right 
in the limiting case of free fall, and since the experimental devia- 
tion from theory is apparently increasing with decreasing oC — the 
trend can also be detected for d^ = 2.5 in. -- we conjecture that 
separation still occurs for contours B and C although these runs re- 
mained turbulent much longer than for contour A. No elementary 
theory can account for this as it involves separation of either or 
both turbulent and laminar unsteady boundary layer on a curved wall. 


It might be possible by a better nozzle design, e.g., adapting the 

mofhnrl nf Thufli'f-pc 11 flnu t'h'rniip-hrmt' 




the run. 


Several runs were made for the nozzle of contour A with 
Z 2 = 13.15 in., 6 = 0.0943. The agreement with theory was essen- 
tially the same as for 6 = 0.0497, and the data are omitted from 
the graph. It is important to note, however, that even for a 
relatively short nozzle the factor 1 - € 2 .ay be significantly less 
than unity. 

CONCLUSIONS 

From the two configurations investigated there is strong evidence 
that unsteady, inviscid, liquid flow from a reservoir, with a free 
surface, should be treated by quasi-steady theory only when varia- 
tions in the cross-sectional area of the reservoir are restricted 
essentially to the neighborhood of the exit station. 

On the other hand it appears that the discharge time of such a 
flow can be determined quite accurately by one-dimensional theory if 
the acceleration of the free surface is neglected. Integration of 
the resulting equation is called the approximate unsteady solution. 

The approximate solution is found not to differ significantly from a 
numerical solution of the exact equation for the conical reservoir, 
whereas for certain ranges of the height ratio it may deviate con- 
siderably from the quasi-steady solution. Agreement between experiment 
and theory for a slender conical reservoir is good. 

Theoretical discharge times for the cylindrical reservoir 
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configuration are about the same for all three theories as long as 
the area ratio is ten or above. For a cylindrical reservoir certain 
precautions are essential if meaningful measurements are to be obtain- 
ed. In the first place to avoid collapse of the free surface it is 
necessary to use a thin, balsa wood disc which rides on the free sur- 
face. Furthermore to eliminate the jet contraction at the exit a short 
transition-section-plus-nozzle must be employed. If the length/dia- 
meter ratio of the nozzle is large compared to unity, and if the Rey- 
nolds number based on the instantaneous outflow velocity is not too 
high, frictional effects may be important. If the ratio is small, 
there appears to be a tendency of the jet to separate from the nozzle 
wall which results in larger experimental discharge times than pre- 
dictable from theory. It ought to be possible by suitable design to 
avoid this separation. Introduction of the nozzle requires a correc- 
tion to the theoretical discharge time which may be large even for 
short nozzles. 

On the question of one-dimensional theory versus the non-steady 
Bernoulli equation, with weighting factors to account for non-uniformi- 
ties, we conclude that the former is the more accurate. In one extreme, 
where experiments' deviate from the one-dimensional theory, the- error 
is due to neglected fluid viscosity which neither theory pretends to 
account for. The other indications are that the discrepancy is due 

'j 

to inadequate design of the nozzle contraction. Thus we do not believe 

that the use of the Kozeny-Stary correction factors can be justified 

for unsteady flows. In fact we suggest that their validity even in 
steady flows should be re-examined. 




30 - 


REFERENCES 


1. Streeter, V.L., editor of "Handbook of Fluid Dynamics," 

1st edition, 1961, McGraw-Hill Book Co. 

2. Roudebush, W.H. and I.I.Pinkel, "Potential-Flow Analysis 

of Unsteady Outflow from a Tank and Its Effect on 
the Dynamics of a Fluid System", NACA TN D-1164, 

June 1962 £ 

3. Prandtl, L. , and O.G.Tietjens , "Fundamentals of Hydro- 

and Aeromechanics , vol. I, p. 120, Dover Publications, 
New York, 1957. 

4. Kozeny, J. , "Hydraulik, " p. 61, Springer-Verlag, Vienna, 

1954 (in German). 

5. Stary, F. , "On the Discharge Time of a Tank with a Horizon- 

tal Bottom Opening," ZAMM, vol. 42, 6, 1962 (in German) 

6. Bird, R.B., W.E. Stewart, and E.N.LightfOot, "Transport 

Phenomena," p. 227, John Wiley & Sons, Inc., New York, 
1963. 

7. Schlichting, H. , "Boundary Layer Theory," 4th edition, 

McGraw-Hill Book Co. , 1960. 

8. Kaufman, Walther, "Fluid Mechanics", McGraw-Hill Book 

Company, Inc., 1963, page 56. 

9. Dwight, H.B., "Tables of Integrals and Other Mathematical 

Data", Third ed. , The MacMillan Co. ,' 1957. 

10. Sabersky, Rolf H. , and Allan J. Acosta, "Fluid Flow", 

The MacMillan Co. , 1964, pp. 65-67. 

11. Thwaites, B. , "On the Design of Contractions for Wind 

Tunnels", ARC Rep. and Memo 2278, March 1946. 




- 31 


APPENDIX 

THE CONICAL RESERVOIR 

The conical reservoir,. Figure 1, was made in two sections. The 
lower part, Figure 8, which consists of three separate pieces bolted 
together, was machined from aluminum round stock, and provided with 
O-ring seals at the horizontal joints. To the upper end of this 
section is attached an extension (not shown) fabricated from 0.060 
inch sheet metal, rolled to the. conical shape. Spot check', cf the 
extension show that the maximum deviation is less than one percent 
across any. diameter at the upper station. The total height of this 
reservoir from the apex of the cone, extended is 3 r inches. 

' To restrain the fluid and to start the flow we use a simple 
‘hinged arm supporting a rubber pad which is pressed against the ori- 
fice lip. The device is held in place by a trigger prior to starting. 
This arrangement works well enough since it does not seem to contri- - 
bute any transients of its own. 

THE CYLINDRICAL RESERVOIR 

For the cylindrical reservoir, Figures 4a and 4b, an empty 
120 lb. commerical grease drum is employed. The sharp-edged orifice 
is obtained by machining a hole in the drum bottom. 

The nozzle configuration uses the same drum with the nozzle of 


■} 
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Figure 9 bolted to the bottom. The nozzle is enlarged by machining 
as we proceed from the one-quarter-inch to the three- inch hole. 


Rubber stoppers are used to cork the supply. 




Reservoirs with A 




Fig. 3 Discharge times for n = 1, 2 
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Fig. 5 Nozzle contours 
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